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A pyrographically decorated gourd, dated to the French Revolution period, has been alleged to contain a 
handkerchief dipped into the blood of the French king Louis XVI ( 1 754- 1 793) after his beheading but recent 
analyses of living males from two Bourbon branches cast doubts on its authenticity. We sequenced the 
complete genome of the DNA contained in the gourd at low coverage (—2.5 X) with coding sequences 
enriched at a higher —7.3 X coverage. We found that the ancestry of the gourd's genome does not seem 
compatible with Louis XVI's known ancestry. From a functional perspective, we did not find an excess of 
alleles contributing to height despite being described as the tallest person in Court. In addition, the eye 
colour prediction supported brown eyes, while Louis XVI had blue eyes. This is the first draft genome 
generated from a person who lived in a recent historical period; however, our results suggest that this sample 
may not correspond to the alleged king. 

The analysis of complete individual genomes is now routinely achieved, providing new data to understand the 
bases for individuality and also to explore the complexity of modern as well as ancient human diversity' In 
the route towards a better characterization of individual genomes and expanding the application of NGS 
techniques into new disciplines, the analysis of genomes from historical periods is the next logical step. 

In 2010, DNA was retrieved from a pyrographically decorated gourd (Fig. 1 and Supplementary Fig. SI) dated 
to the French Revolution period'. According to a text inscribed in the gourd, it contained a handkerchief dipped 
by a witness called MaximUien Bourdaloue into the blood of the French king Louis XVI (1754-1793) after his 
beheading in January 21st, 1793 (Fig. 1). The Y-chromosome haplotype, determined with the AmpFlSTR 
Identifiler PGR amplification Kit is unique in a current European database with over 21,000 individuals, suggest- 
ing that any matching with other individuals would likely support a paternal relationship between them. The 
subsequent analysis of the mummified head" of the king Henri IV (1553-1610)', separated by seven generations 
from Louis XVI, provided a partial Y-chromosome profile that with the exception of one genetic marker, was 
concordant with that found in the gourd's blood. Since one allelic difference is not unexpected in pedigrees of 
several generations, the results still supported that both remains were paternally related". 

However, a later analysis of three male living Bourbons from two different family branches revealed identical 
Y-chromosome haplotypes that were different to those found in the remains of Henri IV and the gourd's 
individual'". To further explore the possibility that the gourd's blood could belong indeed to the king Louis 
XVI and shed further light on these contradictory results, we have retrieved the complete genome of the human 
sample contained within the gourd. The objective was to confirm or deny the authenticity of the sample and to try 
to correlate the genomic data with certain phenotypical traits reported by historical records that might corre- 
spond to the king. 
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Figure 1 | A pyrographically decorated gourd supposed to contain the 
blood of the king Louis XVI (image and copyright by Davide Pettener). 

Results 

Contamination estimates and SNP discovery. We sequenced DNA 
extracted from the inside of the gourd using lUumina technology 
(Supplementary Fig. S2 and Supplementary Methods). Due to the 
inefficiency of the sample, several lanes were used in different 
platforms resulting in a total 2.5-fold genome wide effective 
coverage and 7.3-fold exome coverage (Supplementary Tables Sl- 
S2). About 24% of the sequences were human, 46% corresponded to 
Pseudomonas (a group of bacteria comprising plant pathogens), 27% 
to unknown origin and 1.3% to other organisms, including fungi 
(0.47%) and of course, Cucurbitaceae (0.1%) (Fig. 2) most likely 
derived from the gourd itself The filtered human reads displayed 
the characteristic DNA damage patterns seen at the end of ancient 
sequences, although the signal was lower in comparison to that seen 
in much older samples (Supplementary Fig. S3). 

Considering the highly heterogeneous use of a handkerchief and 
the gourd itself, we anticipated the presence of human contaminants 
in the sample. To estimate this, we made use of the higher coverage 
obtained from the mitochondria (~140X, mtDNA genome cov- 
erage). Different haplo types were reconstructed from the mtDNA 
reads (Supplementary Table S3) with a predominant (73%) Nlbla2 
haplotype, concordant with what had been previously determined by 
polymerase chain reaction (PCR)'; however, three additional haplo- 
types (Hla, Jlc2c2 and K*), not shared with people involved in the 
experimental procedures were present in decreasing ratios (13%, 
9.6% and 1.9%, respectively) adding up to a total upper limit of 
contamination of —24%. 

To obtain nuclear contamination estimates we followed a prev- 
iously developed procedure" that takes advantage that only one X 
chromosome is expected to be present in a male individual (and thus, 
a single allele at each site). The contamination tests are based on the 
analysis of known polymorphic positions (present at 1000 Genomes 
Project Phase 1 data') with low coverage (4X to 10 X) and the com- 
parison of their mismatch rate to adjacent sites (e.g., less than 5 bases 
apart). If the sample is contaminated, the number of mismatches at 
the known polymorphic positions will tend to be significantly higher 
than in adjacent sites. We obtained a 19.4% (CI 95% = 0.179-0.200) 
value of potential nuclear DNA contamination (Supplementary 
Table S4). We also checked the diagnostic haplogroup positions at 
the Y-chromosome finding a similar value of 17% (CI95% = 0.063- 
0.277) of contamination (8 out of 47 reads) (Supplementary Table 
S5). Considering uncertainties associated to sampling and low cov- 
erage, it is likely that contamination is equally prevalent at the nuc- 
lear and at the mitochondrial genome. 



The significant degree of contamination should emerge as hetero- 
geneities in the resulting genotypes. However, we can safely assume 
that in regions with significant coverage, the 17-24% of background 
contamination will consistently emerge as the minor allele. Therefore 
in subsequent analyses we systematically assessed the major alleles. 

To remove contaminants by selecting the major allele, we 
developed a procedure of eliminating SNPs with a substantial skewed 
allele imbalance (Supplementary Tables S6-S7, Supplementary Fig. 
S4). Thus, we finally obtained two sets of single-nucleotide poly- 
morphisms (SNPs): a genome wide set of SNPs with >3X coverage 
(mean coverage = 5X) and a second set of SNPs with >9X coverage 
(mean coverage = 12 X) restricted to the higher coverage exome 
data. In both sets, we used an allele imbalance removal model 
accounting for a contamination ratio of up to 30% (Supplementary 
Methods). While contamination can still remain problematic a lower 
coverages, we have estimated that in the case of 12X coverage the 
remaining contamination, according to a Bernoulli statistical model 
and conditional probabilities, range between 0.25% and 0.39% 
(Supplementary Methods). 

The second set has a higher confidence due to the increased cov- 
erage; nevertheless we have calculated a genotype concordance 
between both datasets of 98%. Obviously, the conservative procedure 
selected had the side effect of removing real heterozygotes from the 
gourd's genome. Subsequent Sanger sequencing of 9 sites confirmed 
the undercalling of heterozygous SNPs (false negative rate of 33.3%) 
(Supplementary Table S8). Combining the two datasets in a non- 
redundant set of SNPs, we characterized a total of 1,208,005 SNPs 
contained in dbSNP, from which 314,549 (26.0%) were reported as 
heterozygous. Only in the exome, we reported 36,928 SNPs from 
which 43.75% are heterozygous. The overall Ti/Tv ratio is 1.91 
(Supplementary Table S6). 

We finally reasoned that our power to make functional interpreta- 
tions from alleles at individual level could still be compromised; 
therefore, we applied a final filtering strategy, based on known link- 
age disequilibrium (LD) patterns in European populations. It is 
expected that contaminants (we known that there are at least three 
within the sample) will generate mixed haplotypes, thus breaking the 
LD blocks on our genome as compared to present-day Europeans. 
We checked if a random sample of haplotypes found in 60 CEU 
(Northern and Western European Ancestry) individuals from the 
1000 Genomes Project pilot phase" were shared with our genome 
(Supplementary Methods). We found 17.898 haplotypes (involving 
—250,000 SNPs) composed by the same linked variants, ranging in 
number from 3 to 85 1 SNPs, and showing an average similarity to the 
described CEU haplotypes over 93.5% (Supplementary Fig. S5). 
Moreover, the haplotypes in the gourd's genome tend to be among 
those found in higher frequencies in modern Europeans (average 
frequency over 60%) (Supplementary Fig. S6). These observations 
suggest that despite the contamination of the sample it is still possible 
to obtain reliable information on the prevalent genome present in the 
gourd. 

Ancestry inference. According to the historical records, the 
genealogical track of the king Louis XVI was highly admixed, with 
a major Central European contribution (Supplementary Table S9). 
His 16 great-great-grandparents were from present-day Germany (N 
= 8), Austria (N = 1), Poland (N = 4), Italy/France -House of Savoy- 
(N = 1), and France (N = 2, one corresponding to Louis, Grand 
Dauphin of France (1661-1711)) (Supplementary Table S9). 

We first determined the gourd's blood paternal ancestry, using 
diagnostic SNPs in the Y-chromosome, as defined by ISOGG 
(2014). Several derived SNPs define a G2a2a haplogroup 
(Supplementary Table S5); again, this result is concordant with that 
previously found with the AmpFlSTR Identifiler kit'. 

An identity-by- descent (IBD) tract sharing analysis showed very 
few IBD tracts with extant Europeans, likely due to gaps in the 
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Figure 2 | Metagenomics analysis and species composition of the gourd's sample. 



genomic coverage and undercalling of heterozygous sites. The 
sequenced individual shared the highest number of fragments with 
a French individual, but shared the longest tract with a Belgium 
individual (Supplementary Table SIO). 

We also performed a principal component analysis (PCA) on the 
>9X SNP exome dataset, using 236 individuals from the 1000 
Genomes Project' for comparison (16,635 SNPs). To minimize dif- 
ferences between our low coverage genome and the reference geno- 
types, we randomly sampled a single haploid allele from each 
individual prior to PCA analysis". The gourd's individual can be 
found close to the CEU distribution, slightly displaced towards the 
Iberian samples and Tuscan samples from Northern Italy (Fig. 3). 
Subsequently, to place the gourd's genome in a more precise 
European context, we performed a genome wide PCA using the 
European populations of POPRES". Because only 3% of our exome 
data intersect with POPRES markers, we utilized all SNPs passing 
our filters genome-wide in this analysis -SNP dataset increased to 
about 100,000. We found the sequenced individual clusters with 
Northern Italian individuals (Fig. 4a) which is not what we could 
expect owing the known ancestry of Louis XVI. This data set is 
obviously more influenced by the contamination of the sample, 
due to the less efficient allele imbalance removal associated to the 
low coverage. To explore if Louis XVI ancestry at the level of his 
great-great grandparents was compatible with the PCA results, we 
have randomly generated 30 composite genomes with the propor- 
tions of his known ancestry and generated a PCA with POPRES 
European populations. This simulation showed that the expected 
position of the king in the PCA is among present day Central 
European populations (Germany and Poland) (Fig. 4b). 

In summary, the results of these analyses do not support the royal 
identity of the sequenced genome. However, given Louis XVI's com- 
plex great-great-grandparental ancestry, the low coverage and the 
aggressive filtering to remove minor alleles from the current dataset, 
we consider it could stiU be possible, although implausible, that the 
gourd's blood could be that of the French king. Alternatively, it could 
be that one of the so-called contaminants detected at mitochondrial 
level may in fact correspond to the king's blood; however, we do not 
have any evidence (e.g., Bourbon Y-chromosome SNPs among our 
reads) that could support this assumption. 

Functional exome assessment and phenotypical traits. We per- 
formed a functional characterization of the variants found in his 



genome, including a catalogue of features unique to this individual. 
A total of 3,850 genetic variants not described in dbSNP were found 
in the exome (Supplementary Table Sll). In the European 1000 
Genomes, this figure is on average 16,059 (min = 15,549, max = 
16,526); the difference can be attributed to the low coverage and the 
contamination removal filtering we applied to the gourd's genome. 
Despite of this, we still found 13 stop codons and 8 non-synonymous 
changes (all in heterozygosity) within genes that are associated with 
known human Mendelian diseases, as well as 7 stop codons and 10 
non- synonymous changes in genes that are not, but occurred as 
homozygote for the alternative allele (Supplementary Tables S12- 
S15). None of these variants could be easily related to any known 
phenotype -including height- in the king. 

We then reviewed different contemporary historical sources, 
including Marie- Antoinette's correspondence and transcribed con- 
versations from testimonies and relatives (Supplementary Methods) 
to characterize physical and behavioural traits of Louis XVI. We 
subsequently contrasted the sequenced genome with several pheno- 
typical features that the king displayed (blue eyes, tall height, obes- 
ity), as well as some possible disorders such as type-2 diabetes. The 
eye colour has an already well established genetic background'''"', 
and tools such as IrisPlex" can yield an accurate prediction of eye 
colour based on the six most-informative single nucleotide poly- 
morphisms (SNPs). 

The remaining phenotypes have complex genetic backgrounds 
that are currently being investigated from genome-wide association 
studies (GWAS). Due to the gaps and reduction of heterozygotes 
associated to the low sequence coverage and the background con- 
tamination, we don't have information for all the published SNPs; 
however, we took advantage of linkage disequilibrium (LD) patterns 
to recapture SNPs not genotyped and also to confirm those present in 
the genomic draft. 

The conservative final set of SNPs after selection and filtering 
(Table 1) did not allow us to screen for type-2 diabetes and obesity 
since virtually no data was available to evaluate the risk susceptibility 
of the individual's genome. Therefore, we continued the analysis on 
two of the traits, eye color and height. 

The king is depicted with blue eyes in all his portraits, despite the 
fact that his parents (Louis, Dauphin of France, and Maria Josepha of 
Saxony) were brown-eyed. In the previous study, we retrieved 
by PGR the critical rsl2913832*C SNP at the HERC2/OCA2 locus. 
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Figure 3 | Principal component analysis (PCA) of the gourd's genome SNP data at >9 x coverage and compared to European individuals of the 1000 
Genomes Project. The analysis was based on 47 randomly selected individuals per population (FINS, CEU, GBR, TSI), except for the Iberians, where all 
six available individuals were included. A total of 32,701 SNPs were included after filtering for heterozygous positions in order to account for the 
background contamination. 
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Figure 4 | (a) PCA of the merged POPRES dataset and the gourd's genome data. For POPRES individuals, abbreviations and colors for each sample 
foUow conventions of a previously published analysis^* in which two letter codes represent population identifiers. The black dot represents the individual 
data from the gourd's genome. A total of 101,107 SNPs (>3X coverage) were included, (b) PCA of the merged POPRES and 30 randomly composite 
genomes generated from the known ancestry of Louis XVI at the great-great grandparents generation. For POPRES individuals, abbreviations and colors 
for each sample follow conventions of a previously published analysis^* in which two letter codes represent population identifiers. Composite individuals 
are represented by uppercase C. A total of 196,526 SNPs were included. 
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SNPs in GWAS catalog 
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After filtering SNPs 


Height 
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58 




Type 2 Diabetes 


1 14 
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finding the nucleotide substitution associated to blue eyes in hetero- 
zygosity. To further explore the gourd's genome phenotype we deter- 
mined the allelic state at five additional loci and predicted the eye 
color with Iris-Plex'^ (Supplementary Table S16). The fact that most 
of the remaining SNPs are associated with brown eyes resulted in a 
low probability for blue (2.4%) or intermediate (8.3%) eye color. 

Complex traits such as height are more difficult to assess from 
genetic data alone. However, we could evaluate 23 high-quality SNPs 
influencing this trait (Table 2). We assumed an additive multiple 
regression linear model for the increase in centimeters that each allele 
adds to the final individual's height (based on the Beta-coefficients 
compiled from literature). To put in context the gain in height 
inferred for the sequenced individuals, we calculated the distribution 
of centimeter gain that each combination of height-increasing alleles 
would confer over an undetermined basal height in Europeans (per- 
forming 100,000 resamplings from CEU allele frequencies under 
Hardy- Weinberg assumption). The allelic combination found in 
the gourd's genome generates an absolute increase of 7.22 cm. 
Given that the average increase in the population is 9.05 cm (SD 
= 1.03), this individual is located in the 3.9 bottom percentile. We 
repeated the analysis using only the genome-wide significant SNPs 
(P < 5 X 10"' and <5 X 10"" in GWAS Catalog), but the individual 
was placed at the 8.37 and 9.78 bottom percentiles, respectively. 
Similar results were obtained when, rather than assuming HW equi- 
librium, resamplings were made from the genotypic frequencies in 
extant CEU samples. In this case, the absolute gain in height for this 
individual is 5.94 cm, which placed the gourd's individual in the 4.71 



bottom percentile (CEU average 6.52 cm; SD = 0.77). The additional 
search for possible mutations in genes associated to cases of gigant- 
ism produced no results. 

Discussion 

We have retrieved a draft genome from the gourd's blood, supposed 
to be that of king Louis XVI, and we have conducted several genomic 
analyses to characterize this individual. First, we have generated a 
population analysis with different current European datasets as ref- 
erence populations. These analyses show that the position of the 
gourd's individual in the context of modern Europeans is somewhat 
complex - we see IBD sharing with contemporary French indivi- 
duals, but also ancestry that is somewhere between contemporary 
CEU samples and Northern Italian samples. Although it can be 
argued that Northern Italian populations have clear ancestral ties 
with Central Europe, this discernible component cannot be easily 
reconciled with the known Louis XVI family history, given that just 
one in sixteen great-great-grandparents - Victor Amadeus II, Duque 
of Savoy (1666-1732)- has some possible Northern Italian ancestry. 
Additionally, it has to be remembered that the gourd itself has been in 
Italy for more than a hundred years, but we have no evidence that the 
contamination found in the sample could be traced back to this 
period. 

Secondly, we have provided a functional interpretation of the 
gourd's genome, both at the level of disease susceptibilities and also 
physical traits. Due to the limitations associated to the low coverage, 
we have focused on the study of height and eye colour. 



Table 2 | Height associated SNPs ascertained for the phenotype inference analysis from the gourd's genome. The allelic load column 
represents the number of height increase alleles for that SNP that our individual carries. EAF ("Effect Allele Frequency") is the allele that 
exerts an increase in cm in height. The allelic load can beO, 1 or 2 depending if the individual carries none, one copy or two copies of the EAF 
allele, respectively. Mean refers to mean in height of the individuals analyzed and SD is the standard deviation 

gourd's DNA 

chr. region position pubmed ID SNP EAF Beta/Allele (cm) pvalue mean SD allelic load 
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Contemporary witnesses stated that the king was the tallest person 
at Court, which was arguably a representation of the best fed people 
in the country. Some estimates, including one deduced from the 
length of the coronation cloak in 1775 (measuring 162 cm), placed 
Louis XVI as an extremely tall person, around 185-190 cm in height 
or even more"*. Our estimates suggest the genetic background of the 
gourd's genome could have predisposed to a height increase of 
~7 cm, but this increase is in the lowest range of height increase 
expected for random present-day individuals of European ancestry. 
Thus, this individual's height should be only slightly over the average 
male height in the country that in 1750, as deduced from military 
French records, was of — 167 centimeters'". Despite that some SNPs 
influencing this trait are likely still unknown, no gene-by-envir- 
onmental effects have been described so far for height GWAS 
SNPs. Additionally, height SNPs display a high replicability in other 
human populations, such as Africans and Asians^". Therefore, 
although it cannot be discarded that the perception of the king's 
height would be somehow exaggerated by adulation-driven cour- 
tesan attitudes, or that other factors (e.g., rare variants not captured 
by current GWAS) were partially responsible for the king's height, 
our observation here suggests that the genome is inconsistent with 
that of a very tall person. 

Despite having in heterozygosis a major-effect SNP for blue eye 
colour, four other SNPs involved in this phenotype display only the 
allele associated to brown eyes, to the point that the estimated pre- 
diction probability is 0.892 for brown eyes. Although this could be 
considered an upper limit due to the underrepresentation of hetero- 
zygous positions in four SNPs not located in the exome (currently 
between 1-4X of sequence coverage), the prediction for blue eyes 
could only be increased up to 0.195 if these SNPs were in fact het- 
erozygous. Thus again, within the limit of the quality of our data, the 
sequenced genome appears to be inconsistent with that of a blue- 
eyed person. 

Finally, we have generated a catalogue of private genetic features, 
including stop codons and non-synonymous substitutions. If 
Bourbon genomes contemporaneous to Louis XVI can be retrieved 
in the future, it would be relevant to the authenticity question to 
determine if the private mutations found in the gourd's genome 
are present in members of the French royal family. 

Although we cannot totally discard that the gourd's sample 
belongs to Louis XVI on our genomic data alone, several lines of 
evidence, including the ancestry analysis and the functional inter- 
pretation of the genome fail to provide definitive support for the 
attribution of this specimen to the beheaded French king. Higher 
quality genomes of this sample as well as that of other possible 
genomes from the Bourbon family will be needed to make a more 
definitive claim. 

Second-generation sequencing techniques will allow the genera- 
tion of massive genomic data in forensic historical cases like those 
associated to relicts from the French royal family. These data grant 
more accurate inferences on ancestry and phenotype than those 
obtained by traditional forensic methods, although some of the dif- 
ficulties of identifying personal samples remain. With this approach 
we have obtained and analysed here the first draft genome from a 
person of a recent historical period. Irrespective of the historical 
significance of the individual analysed, this kind of genomic data 
can provide a new source of information on disease susceptibility, 
demography, mobility patterns and ancestry in historical periods. 

Methods 

DNA extraction, sequencing and mapping. About 60 mg of dried blood sample 
obtained from the inside of the gourd was extracted in a dedicated ancient DNA 
laboratory in Barcelona, following experimental procedures described elsewhere'^\ 

We generated in several phases a total of 1,545,311,190 reads with Illumina tech- 
nology (Supplementary Methods). After removing residual primer sequences -a 
problem expected from the sequencing of short ancient DNA fragments^, 214,531,257 
reads were mapped to the human reference genome (hgl9). PGR duplicates and low 
quality reads were also removed leaving a final set of 51,507,079 uniquely mapped 



sequences which correspond to a fmal genome-wide coverage of — 2.50X. We also 
used lUumina technologies to sequence the exome (64 Mbps) captured using 
NimbleGen in-solution capturing technology. A total of 277,762,784 sequences were 
mapped, from which 9,565,693 remained after removing duphcates, thus achieving a 
final global coverage of "-7.32X. 

SNP assessment from Linkage disequilibrium patterns. To recapture non- available 
SNPs and to confirm genotyped SNPs, we first searched for those SNPs present in the 
NHGRI GWAS Gatalog"*^ through a custom-designed pipeline. As a primary input, 
we intersected the HapMap SNPs (CEU population) with the SNP positions analysed 
in the gourd's genome. The rationale behind is that any SNP in HapMap would have a 
frequency high enough to be detected in a GWAS study. This resulted in an input 
dataset of 2,031,376 SNPs that were further analysed to ascertain LD patterns among 
them. 100 kbp regions around each ofthese—2 million SNPs were defined and all the 
SNPs with a LD threshold of R** — 1 were considered. We then searched those 
recaptured SNPs associated to our selected traits at the GWAS Catalog. Afterwards, 
based on sequence data for the CEU population from the 1,000 Genomes Project', 
pairwise haplotypes with recaptured and genotyped SNPs were built using PLINK- 
1.07'^''. With the pairwise haplotypes we could infer if a particular non-genotyped risk 
allele is segregating in LD with any of the known positions. We subsequently filtered 
the obtained SNP dataset and considered only those SNPs recaptured by at least three 
neighboring SNPs, with no incongruities (Supplementary Methods). 

Furthermore, the pipeline was used not only to recapture unknown SNPs, but also 
to verify that the SNPs for which we already had a genotype were correct. We 
considered as valid those SNPs that were in a congruent haplotype with at least six 
other SNPs (Supplementary Table S17). Finally, we checked for LD between the 
recaptured SNPs, to identify putative duplicated tagged regions. We considered that 
two SNPs constitute a single signal if the LD between them had an R^ > 0.8, keeping 
the one with the lowest p-value. 
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